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A, recently presented, general procedure to find static and axially symmetric, interior solutions to 
the Einstein equations, is extended to the stationary case, and applied to find an interior solution 
for the Kerr metric. The solution, which is generated by an anisotropic fluid, verifies the energy 
conditions for a wide range of values of the parameters, and matches smoothly to the Kerr solution, 
thereby representing a globally regular model describing a non spherical and rotating source of 


gravitational field. 


In the spherically symmetric limit, our model converges to the well known 


incompressible perfect fluid solution. The key stone of our approach is based on an ansatz allowing 
to define the interior metric in terms of the exterior metric functions evaluated at the boundary 
source. The physical variables of the energy-momentum tensor are calculated explicitly, as well as 
the geometry of the source in terms of the relativistic multipole moments. 


PACS numbers: 04.20.Cv, 04.20.Dw, 97.60.Lf, 04.80.Cc 


I. INTRODUCTION 


Since the discovery of the Kerr metric there have 
been many attempts to find a physically meaningful mat- 
ter distribution that could serve as its source (see 
and references therein). However, no satisfactory solu- 
tion has yet been found to this problem. 


It is the purpose of this work, to provide an interior so- 
lution to Einstein equations, satisfying all the usual phys- 
ical conditions, and smoothly matched on the boundary 
surface of the fluid distribution, to the Kerr metric. With 
this aim we shall generalize a procedure, recently pro- 
posed to find sources of the Weyl metrics IER to the 
stationary case. 


As a particular example we shall find a source for the 
Kerr metric which consists in an anisotropic fluid, satis- 
fying the Darmois matching conditions on the boundary 
surface of the matter distribution, thereby excluding the 
presence of thin shells, and exhibiting a well behaviour of 
all physical variables, for a wide range of values of the pa- 
rameters of the solution. These include values which are 
commonly assumed in realistic models of rotating neu- 
tron stars and white dwarfs. 
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II. THE GLOBAL MODEL OF A 
SELF-GRAVITATING STATIONARY SOURCE 


A. The exterior metric 


The line element for a vacuum stationary and axi- 
ally symmetric space-time, in Weyl canonical coordi- 
nates may be written as : 


ds% = e°” (dt wd)? e 20 +20 (do*+dz7)+e 2% p*d¢*, 
(1) 
where y = W(p,z) , T = T(p,z) and w = w(p,z) are 
functions of their arguments. 
For vacuum space-times, Einstein’s field equations im- 
ply for the metric functions 


Foot fot f) f- A 


f(W po + p W, p + Wizz) + 2f pwp + 2f,2w,2=0, (3) 
with f =e?” and 


I ae E 
T = 1 Cie = fa _ 7’ oF _ w?,) 
T e I 
I = z2 “Tolz a 5? |Pw piles (4) 
Notice that (2),(B) are precisely the integrability condi- 
tion of (4), that is: given any w and w satisfying (2),(), 
a function T satisfying always exists. 
We can write the line element above, in Erez-Rosen 
E9), or standard Schwarzschild-type coordinates {r, y = 
cos 0} or in spheroidal prolate coordinates {x = —“, y} 


ba. M 


pe=r(r—-2M)(1—-y’?), z=(r—M)y, (5) 


where M is a constant which will be identified later. 
In terms of the above coordinates the line element (I) 


| 


may be writen as: 


ds% = —e?t y) (dt — wd¢)? + eo 2P+ 2 ry) -T gp? 4 eo 2lb—¥°]4+21P Oy) —P'] 2 gg? 4. e—21¥—- "lp? sin? 0do?, (6) 


where Y° and [° are the metric functions correspond- 
ing to the Schwarzschild solution, namely, 


1 (==) 1 a 
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Pe oa r(r — 2M) : 
(7) 


where the parameter M is easily identified as the 
Schwarzschild mass. 


B. The interior metric 


We shall now assume for the interior axially symmetric 
line element: 


2g—2a ae 
dr? ait E29 24,2 Gg 


ds? = -Z (r)? (dt — Odd)? + 


+ e747? sin? Odd’, (8) 


ĝ = g(r, 0) — g°(r,0), (9) 


where Q = Q(r,0), and a(r) and g°(r,0) are functions 
that, on the boundary surface, equal the metric func- 
tions corresponding to the Schwarzschild solution (7), i.e. 
af (rz) = Y$ and g°(rg) =T$. Also, A(r) = 1 — pr? and 


1 
AZ = Arx) — zV AC), where p is an arbitrary con- 


stant and the boundary surface of the source is defined 
by r = rg = const. 

The case w = 0, ĝ = â = 0, corresponds to a spher- 
ically symmetric distribution, more specifically, to the 
well known incompressible (homogeneous energy density) 
perfect fluid sphere, and hence the matching of (8) with 


the Schwarzschild solution implies p = —,-. The simple 
Fo 


condition w = 0 recovers, of course, the static case. 

It should be noticed that for simplicity we consider here 
only matching surfaces of the form r = ry = const, of 
course more general surfaces with axial symmetry could 
however be considered as well. 


C. The matching conditions 


We shall now turn to the matching (Darmois) condi- 
tions . Thus the continuity of the first and the second 
fundamental form across the boundary surface implies 


m 


the continuity of the metric functions and the continuity 
of the first derivatives 0,911, 0-900, Oras, producing: 


azs =Ņġs, a:=%y, ge=Tn, gs=T5, 
a=, (a)s = (Y°)s, 

gs =TS, (9*)5 = T°)s, 

Qs =ws, W=ws, (10) 


where prime denotes partial derivative with respect to r 
and subscript © indicates that the quantity is evaluated 
on the boundary surface. It is important to keep in mind 
that we are using global coordinates {r, 0} on both sides 
of the boundary. 

Thus, our line element (8) matches smoothly with any 
stationary exterior (6), provided conditions (IQ) are sat- 
isfied. 

In the particular case when we want to match our in- 
terior with the Schwarzschild exterior (the static limit), 
then Y = ~* and I = I®, and the source is a perfect fluid 
ifa=g=0N=0. 

We shall now see how the field equations constrain fur- 
ther our possible interiors. 


D. The field equations and constraints 


Let us first analyse the well known case when the in- 
terior is spherically symmetric, then â = g = 0, and the 
physical variables are obtained from the field equations 
for a perfect fluid, the result is well known and reads (in 
relativistic units) 


-=u = — 


T =T =T=P = 


2m(r) 


with A = 1— =1—pr=1- , where u 


and P denote the energy density and the isotropic pres- 
sure respectively, and for the mass function m(r) we have 


mr) = -ar | r°Tldr, (12) 
0 


implying 


>) pre 
M=m(ry) = -4r | rid = =. (13) 
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This model, which describes the well known incom- 
pressible perfect fluid sphere, is further restricted by the 
requirement that the pressure be regular and positive 
everywhere within the fluid distribution, which implies 


r 9 
= > =. As it is evident from (11) the pressure van- 


ishes at the boundary surface. 
Finally, if we impose the strong energy condition P < 
H, we should further restrict our model with the condition 
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We shall now proceed to consider the general, non— 
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with respect to y = cos 0, and 
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spherical case. Thus, for our line element we have 
the following non vanishing components of the energy— 
momentum tensor: 


-TÌ = «(8x + Pz — E +365, +600), 
Ti = «(8%P — fre —6J_), 
Tz = «(8"P + Prs +ÔJ-), 
T = k (8TP — pzz + 601+ 65,), 
TÈ = —«06l, (14) 
(1 — A) ` R Q 9’ 
rVA(3VAx — VA) eee sa) re sin? 0] ’ (15) 
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From the expressions above, using (14)-(16) and intro- 
ducing the dimensionless parameter s = r/ry, we can 
now obtain the explicit expressions for the physical vari- 
ables, (with A = 1 — (2s?) /T, Ay = A(s = 1)): 


kt {6 â 9.06 9.5 Ge. 2 cos 0 
-TP = £22 AG Pr + (ôs) (âs) a ee a t â o aeae 
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Obviously for any specific (non-spherical) model we need 
to provide explicit forms for â, g and 2, however even at 
this level of generality we can assure that the junction 
conditions (10) imply (P,. = grrT})s = 0. 

We shall first proceed to prove the above statement, 
and then we shall provide a general procedure to choose 
â and g producing physically meaningful models. 

It is always possible to choose the metric functions â 
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where dy = ps- y$, Îs = Ts-T$, and dot also denotes 
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(20) 
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(21) 
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and g such that, once the junction conditions (IQ) are 
satisfied, the angular derivatives of such functions are 
continuous, i.e.: (G.9)s = (w.9)s and (G,0)5 = (T o)s, as 
well as (Q os = (we)s . 


2M 
in (4), we obtain for 


xy 
Dex + 6J— on the boundary surface: 


Then, using Ay = = 


+ 2M be + re(rs —2M)b2 — (rs — M)Îh+ 
22 
= (22) 
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partial derivative with respect to 6. 


Taking into account the Einstein’s vacuum equations 
(Gag = 0), we find the following relation for the deriva- 
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—ret¥ w? 
4(r — 2M) sin? 6 |r? 


implying the vanishing of (T})s. 
In a similar way it can be shown that T? vanishes on 
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et! = Ar(r — 2M) yp! 


tives of the metric function w, from the Einstein tensor 
component G41 = 0: 


the boundary surface, if we take into account that the 
Einstein vacuum equation G12 = 0 produces 


sin? 0 


and from (I5), we have the following expression for T? 
on the boundary 


ty =. ait ay ay COSO 
(Tr)s = 2(e)svs -Ts Sing F 
Try — M x 2M 
a SS 
ry (ry = aM) 0)s T ry (ry — 2M) (0) k 
"o (25) 


2rs(rs — 2M) sin? 6 


This last condition (T?)» = 0, which follows from the 
Darmois conditions, and therefore is necessary, in order 
to avoid the presence of shells on the boundary surface, 
can be obtained at once from a simple inspection of the 
equation (22) in [29]. 


E. The ansatz for the metric functions 


We shall provide a general procedure to choose â, g 
and Q producing physically meaningful models. With 
this aim, besides the fulfillment of the junction conditions 
(10), we shall require that all physical variables be regular 
within the fluid distribution and the energy density to be 
positive. 

To ensure the fulfillment of the junction conditions 
(10), we may write without loss of generality, 


â = dy(0) + by(0)(r — rs) + A(r,0)(r — ry)? 
g = Îs(0) + Î5(0)(r — rs) + E(r, 0)(r — rs)’, 
Q = ws(0) + wy(0)(r — rs) + I(r, 0)(r — rs)?,(26) 


where A(r, 0), E(r, 6) and I(r, 0) are so far two arbitrary 
functions of their arguments. 

On the other hand, to guarantee a good behaviour of 
the physical variables at the center of the distribution we 


, (24) 


shall demand: 


âh = Â 90 = Â 900 A’ 9 a 560 0 
ĝo = 9,00 = 9,000 = G90 = G'o90 = 9 
ĝo = ĝ'oo = O, (27) 
Q = 0% = 0% = 07" =0 
Q= =0 
Ope Oia OY = OF eG, (28) 


where have been used, and the subscript 0 indi- 
cates that the quantity is evaluated at the center of the 
distribution. 

Using the conditions above in (26) we may write for A, 
= and II 


A(r,0) = Ao(8) + Ab(0)r + F(r, 8) 
E(r,0) = Eol) + So(0)r + 2h (0)r? + G(r, 8) 
I(r, 0) To(0) + I6 (0)r + 11g (0)r? + 
+ Ij’ (@)r? + H(r, 0) (29) 
with F(0,0) = F’(0,0) = 0 and G(0,0) = G’(0,0) = 


G” (0,0) = 0, as well as H(0,0) = H'(0,0) = H” (0,0) = 
H” (0,0) = 0. Then we can finally write for â and g 

lease notice a that there is a misprint in Eq.(43) in 
T there it should read I instead of w in the two terms 
within the round brackets for the metric function g(r, 0)) 


a(r,0) = r° (5 + 2) +r’ (4 — 2) 
+ (r—ry)?F(r, 6) (30) 


+ (r—ry)?H(r, 6). (32) 


The metric functions obtained so far, satisfy the junc- 
tion conditions and produce physical variables which 
are regular within the fluid distribution. Furthermore the 
vanishing of g on the axis of symmetry, as required by 
the regularity conditions, necessary to ensure elementary 
flatness in the vicinity of the axis of symmetry, and in 
particular at the center (see 23}, (24h, Ba). is assured by 
the fact that Îs and RA vanish on the axis of symmetry. 
Furthermore, the good behaviour of the function 2 on the 
symmetry axis is fulfilled since wy and wi, vanish when 
y = +1. Finally, let us note that the energy-momentum 
tensor components (14), do not diverge on the sym- 
metry axis because the first derivative with respect to 
the angular variable 0 of both Q’ and Q” vanishes on the 
symmetry axis since not only wy and w vanish there, 
but their first derivatives with respect to 0 vanish as well. 

So far we have presented the general procedure to build 
up sources for any stationary metric, in what follows, 
we shall illustrate the method with the example of Kerr 
metric. 


II. PARTICULAR SOLUTION 
A. <A source for the exterior Kerr’s solution 


The Kerr metric in Weyl coordinates is given by the 
following metric functions 


(rı re (1 — 97) — 4M — 97) + rr — r2)? 


(rı +r2 +2MP(1- 3}?) +3? (ri =r)? 
(33) 


f= 


op _ (rı +r (1 -— 3J?) — 4M? (1 -— 9?) +3? (ri — ro)? 
e E Se O 
4rıra (1 — j?) 
(34) 


j(2M +71 +r2)(4M?(1 — j?) — (rı — r2)?) 

(1 PU =F) — MP FP — 

(35) 
where j = = =a/M denotes the dimensionless param- 
eter representing the angular momentum of the source, 
and is related to the rotation parameter a of Kerr in 
its well-known Boyer-Lindquist representation. Also, 
r12 = y p? + (z + My1- j?)?, which in the Erez-Rosen 


coordinates becomes 


2 
ra = |(r - M) + My(r - M)VTZ P| —M?;?(1-y). 
(36) 
As is well known, the relativistic multipole moments 
(RMM) [2o80], of the Kerr solution, are easily given 
in terms of the rotation parameter j (as well as a). In 
fact, using the FHP method in order to calculate 
the RMM, these are expressed in terms of the expansion 
coefficients (mọ) of the Ernst potential on the axis of 
symmetry, but the Kerr solution is the only one verifying 
that each RMM (Mpx) at any order k is just equal to 
the corresponding coefficient mg for such order, which 
implies, 


Therefore, the massive RMM (even orders) and the ro- 
tational RMM (odd orders) can be expressed as follows 


Ma, = ae Mai+1 = Maly Mre jai, 
(38) 

In particular let us remind that a first conclusion de- 
rived from these RMM (8), is that the rotation of the 
object leads to a negative quadrupole massive moment, 


M 
q= Te = —j*, i.e. all the possible sources of Kerr solu- 
/ 
tion are oblate. 


Let us next, consider the line element with metric 
functions given by (80) and (81) 
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a(r,0) = ýs? (3 — 28) + rsýhs?(s — 1), 
Pys(4 — 3s) + rgs? (s —1), 
Q(r,0) = wyst(5 — 4s) +rswhs (s — 1). (39) 
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with s=r/ry € [0,1]. 
For the Kerr solution we may write: 


ds = by — y$ = 


Pelr Š= 


eE oe 


wy = Mj 


=—r(r-2)+(1-y’) -j?. (42) 

A straightforward calculation, using ([7)—(21) allows 
us to find the explicit expressions for the physical vari- 
ables, these are displayed in figures (3)-(7). However, 
before entering into a detailed discussion of these figures, 
we shall carry out some calculations with the purpose 
of providing some information about the “shape” of the 
source. In particular we shall see how it is related with 
the rotation parameter 7 of the source. 

For doing so, using our interior metric, we shall calcu- 
late the proper length l, of the object along the axis z, 
and the proper equatorial radius lp: 


i F ge i r Aa a 
z= TSBs = ay a a 
0 VA do VA 
(43) 


where p,z are the cylindrical coordinates associated to 
the Erez-Rosen coordinates. 

Obviously in the spherical case (â 
lengths are identical: 


p= f° er =Ty V5 arcsin = (44) 
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where the fact that p = — > has been taken into account 
Tr 


and where l$, |; denote the lengths corresponding to the 
spherical case. 

It would be convenient to introduce here the concept 
of ellipticity (e), which in terms of l, and lp, is defined as 
e=1- t, The two extreme values of this parameter are 
e=0, which corresponds to a spherical object, and e = 1 
for the limiting case when the source is represented by a 
disk. In between of these two extremes we have e > 0 for 
a prolate source and e < 0 for an oblate one. 

In the general (non-spherical case) we must compare 
function e~@4=) with e94¥=9)-44=9) since g(y = +1) 
vanishes along the axis. It can be seen that the sign of 
both functions is positive, and their relative magnitudes 


M sae 
(7-1)? — y? N+ rP rz (27? — 1) 
nd rr —2) PFG) }, 
(N+ rPr3)\(2 +r? +73) 

(J? — 1) [G2 —1)(N + rPr}) + 2(N + 94)’ 


N + rr (27? — 1) } 
— P)(rp + rz + 2) 


(40) 
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verify the inequality e7451) > e9(¥=9)-4(y=9) for all val- 
ues of s in the range s € [0,1], no matter the sign of the 
parameter j, leading to the well known result that the 
rotation of the object always generates an oblate source 
(q < 0) since l; < lp. In the figure 1, one example is 
shown. 
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FIG. 1: Functions exp(g(s,y = 0) — a(s,y = 0)) and 
exp(—a(s,y = 1)) for values of the rotation parameter j = 
+0.1 and T = 2.7. 


Figure 2 shows the ellipticity e of the source as a func- 
tion of the rotation parameter j, for different values of 
the parameter T. As can be seen, the relation between 
e and j for any value of r shows that the greater is j, 
the greater is e, and therefore the shape of the source is 
more oblate. Of course, for j = 0 (static case) we recover 
the sphericity (e = 0). It is also observed from the fig- 
ure 2 that the deformation of the source with respect to 
the spherical case, for any fixed value of j, is smaller for 
larger values of 7 (less compact object). 

Let us now turn back to the physical variables of our 
model. Figure 3 exhibits the behaviour of the radial pres- 
sure Pr. = grrT} for different values of j. 

In it, we observe the variation of the radial pressure 
with respect to the spherically symmetric case (j = 0). 
This variation is smaller for angle values close to the 
equator, as it is apparent for y = 0.3. Notice that the ra- 
dial pressure is positive, with negative pressure gradient, 
and vanishes on the boundary surface. 

Figures 4 and 5 depict the behaviour of different en- 
ergy momentum components, for a specific choice of the 
parameters q and 7T. 

Figure 6, shows the verification of the strong energy 
condition (—Tj) — T} > 0 (also for a specific choice of 
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FIG. 2: Relation between the ellipticity of the source e and its 
rotation parameter j for different values of T. 
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FIG. 3: Four profiles of r&P,r = r&grrTt, as function of s, 
for y = 0.3 (graphic a), and y = 1 (graphic b) with T = 2.7, 
and different values of j. 


the parameters j and 7). 
Figure 7, shows the component T8. 


IV. DISCUSSION 


By extending the general procedure developed in 
to the stationary case, we have been able to build up a 
physically meaningful source for the Kerr metric, satisfy- 
ing the matching conditions on the boundary surface of 
the matter distribution. In spite of the fact that_a per- 
fect fluid source for the Kerr metric might exist Bd, our 
source is necessarily anisotropic in the pressure. 

Particular attention deserves the presence of a non- 
vanishing Të component of the energy-momentum ten- 
sor. Indeed, defining as usual an energy-momentum flux 
vector as: F” = —V*T,,, (where V” denotes the four 
velocity of the fluid), it appears that, in the equatorial 
plane of our system, energy flows round in circles around 
the symmetry axis. This result is a reminiscence of an 


FIG. 4: —r&T§ (graphic a), and r£T (graphic b), as func- 
tions of y= cos@ and s, with j = 0.1 and T = 2.7. 


FIG. 5: (a) r2,.T3,and (b) rT? as functions of y and s. 


effect appearing in stationary Einstein—Maxwell systems. 
Indeed, in all stationary Einstein-Maxwell systems, there 
is a non vanishing component of the Poynting vector de- 
scribing a similar phenomenon (33), Bal (of electromag- 
netic nature, in this latter case). Thus, the appearance 
of such a component, seems to be a distinct physical 
property of rotating fluids, which has been overlooked 
in previous studies of these sources. 

We have carried out a systematic search for the range 
of values of 7 and j for which our models exhibit accept- 
able physical properties. We have focused on the fulfill- 
ment of Positive Energy Density (P.E.D.) (~T? > 0), 
Strong Energy Condition (S.E.C.) ((—T?) — T? > 0) and 
Positive Radial Pressure (P.R.P.) (Prr = g1iT} > 0).The 
results of this search are shown in Table I. It appears ev- 
ident that physically meaningful sources exist for a wide 
range of values of the parameters. 

On the other hand, the chosen range of values of the pa- 
rameters, incorporates values considered in the existing 
literature, to describe realistic models of rotating neu- 
tron stars and white dwarfs. Let us elaborate on this 
issue with some detail. 

The rotation of the source is determined by the param- 
eters a, M of the exterior solution. Indeed, the rotation 
parameter j = J/M? = a/M stands for the dimensionless 


FIG. 6: r3[(—T)) — TH] as function of y and s, with j = 0.1 
and T = 2.7. 


FIG. 7: r3.T3 as function of y and s, with j = 0.1 and T = 2.7. 
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